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Abstract 

A comparative study of various parametrizations of the dark energy equation of state 
w{z) and its variation with the redshift is made. Astrophysical constraints from 'integrated 
tracking' are laid down to test the physical viability and cosmological compatibility of these 
parametrizations. A critical evaluation of the 4-index Hannestad parametrization is done. 
It is found to be simple and transparent in probing the evolution of the dark energy during 
the expansion history of the universe but it does not satisfy the BBN constraint on the dark 
energy density parameter for the best fit values; however it looks more compatible with the 
observations in the extended parameter space with w << —1. 



1 INTRODUCTION 



The recently discovered 16 Type la supernovae with the Hubble Telescope (Riess et al,2004) 
provide a distinct scenario and conclusive evidence of the decelerating universe in the past 
{z > 1) evolving into the present day accelerating universe. Thus the existence of dark en- 
ergy, which accelerates the cosmic expansion , has been firmly established and the magnitude 
of its energy density has been precisely determined (Wang and Tegmark,2004). The focus 
is now on the evolution of dark energy and its equation of state with the cosmic expansion. 
The simplest and the most natural candidate for dark energy is the cosmological constant A 
with a constant energy density p\ and a fixed equation of state parameter w — —1. But the 
recent analysis of the SNe data (Alam et al 2004;Hutcrcr ct al 2004) indicates that the time 
varying dark energy gives a better fit to observational data than a cosmological constant. 

There are two contrasting prevalent views about the evolution of dark energy at present. 
Riess(2004) and Jassal ct al(2004) have argued that the current SNe la observational data is 
inconsistent with the rapid evolution of w{z). On the other hand, Bassett et al(2004) contend 
that the inconsistency arises on account of the inadequacy of 2-index parametrizations. 
They claim that the rapid variation in w{z) in fact provides a better fit to the 'gold set' 
SNe observations, even after including CMB and large scale structure data (Corasaniti et 
al,2004). According to our analysis, the rapidity of variation |^|, apart from other factors, 
depends on the absolute value of \w\. For quintessence models — 1 < w < — |, whereas w 
has no lower bound for phantom models. Hence the equation of state of dark energy varies 
more rapidly for phantom models with large 1^1 than for quintessence models, other factors 
remaining unchanged. The fact that the observational data (Caldwell ct al,2002) including 
SNe and galaxy clustering, shows a bias towards phantom models might explain why rapid 
variation in w{z) provides a better fit to observational data. 

It is difficult to know the exact functional form of the equaion of state w{z) obser- 
vationally as such different parametric forms (Huterer and Turner ,2000; Weller and Al- 
brccht,2001;Polarski et al,2001;Linder,2003;Padmanabhan ct al,2004; Corasaniti et al,2004; 
Alam et al,2003 ) and non-parametric forms (Bassett et al,2002; Corasaniti and Copeland,2002; 
Corasaniti et al,2004; Bassett et al,2004) of u'(2;)have been used to simulate the behavior 
of evolving equation of state of dark energy and comparison made with the SNe observa- 
tions. The SNe observations essentially measure the luminosity distance (/^(z) which when 
compared with the theoretical parametric values, yields the best fit values of the parameters 
(For details see Padmanabhan,2003). Odintsov and Nojiri(2004) have discussed the thermo- 
dynamics of the evolution of dark energy and have stressed that the phantom stage might 
be a transient period in cosmic evolution. 

For making a comparative study of various parametrizations and their viability, we need 
some criterion. In our previous papers (Johri, 2001, 2002), we introduced the concept of 'in- 
tegrated tracking' and outlined certain astrophysical constraints to be satisfied by the dark 
energy fields for realistic tracking. For example, the dark energy density parameter at nu- 
cleosynthesis epoch {z — 10^°) must satisfy the constraint 

(A) {nx)BBN < 0.14 

The latest analysis (Olive et al,2004) constrains 
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(A) inx)BBN < 0.21 

Again at the point of transition from deceleration to acceleration phase of expansion, a 
viable parametrization must satisfy the constraint (corresponding to g = 0). 

In this paper, we have discussed the merits and demerits of four dark energy parametriza- 
tions, first three involving 2-parameters and the fourth (Hannestad, 2004) involving 4- 
parameters. More parameters mean more degrees of freedom for adaptability to observa- 
tions, at the same time more degeneracies in the determination of parameters. Hannestad 
parametrization is critically examined for its effectiveness and limitations to represent the 
variation of w{z) in compatibility with the observations. 

The outline of the paper is as follows In section 2, we have discussed the expansion history 
of the universe. Assuming a spatially flat universe, the field equations give the Hubble expan- 
sion and deceleration parameter in terms of function f{z) which involves integral of the vary- 
ing equation of state of dark energy. In section 3, we have discussed 2-index parametrization 
models by Huterer and Turner(2002),Weller and Albrecht(2002), Linder(2003) and Padman- 
abhan et al (2004). In section 4, we have critically examined Hannestad parametrization 
model which involves 4 parametrs imparting four degrees of freedom to choose them. In 
section 5, we use interpolation technique to study the behavior of the dark energy under the 
assumption of slowly varying equation of state. In section 6, we conclude with some remarks 
on parametrization methods. 

2 EXPANSION HISTORY OF THE UNIVERSE 

Assuming a spatially flat (k=0) friedmann universe, the fleld equations are 

= H',[ni,{l + zf + Q'^^f{z)] (1) 

and 

^ = wiz)nx (2) 

where if = ^ is the Hubble constant, q = — is the deceleration parameter, w{z) = ^ 
is the equation of state of dark energy and fix = 1 ~ is the cosmic dark energy density 
parameter with fix — flxfiz), 

f{z) = e^v{^[^^dz') (3) 
The luminosity distance is given by 

We can test any parametrization by taking the best fit values and applying constraints 
(A) and (B) to check its compatibilty with the cosmological observations. The transition 
redshift may also be derived from the best flt values. 
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3 Two Index Parametrizations 



I. Linear Redshift Parametrization of w{z). 

(Huterer and Turner,2001; Weller and Albrecht,2002) 

w{z) = Wo + w'z, (5) 
Inserting for w{z) in Eq.(3), we get 

Qx = ^xi^ + X exp[3w'z] (6) 

It has been used (Riess et al, 2004) for probing SNe observations at 2; < 1. It diverges for 
large z. The best fit values for SNe 'gold set'(Dicus and Repko,2004) , Wo = — 1.4, w' = 
1.67, = 0.3, zt — 0.39 give i^xizr) — 0.63 whereas constraint (B) demands Qx — 0.445. 
This parametrization gives +40% deviation from the acceptable value of the dark energy 
density. 

II. Polarski and Linder Parametrization (Polarski,2001;Linder,2003) 

w{z) = Wo + Wi—— (7) 
J- ~r Z 

On differentiation, ^ = It indicates rapid variation of w{z) at z=0 which goes on 

decreasing with increasing z. For large z, w{z) = Wo + Wi. 
Eq.(3) gives 

Qx = 0^(1 + ^)'^'+"'°+"'^) xexp[-^^] (8) 

1 + z 

Best fit values for SNe 'gold set' (Dicus and Repko,2004;Gong,2004) Wo = -1.6,^1 = 
3.3, fllf = 0.3, Zt = 0.35 give fij = 0.626 whereas contraint (B) demands fl^ = 0.46; 
as such this parametrization gives +36% deviation from the acceptable value of dark energy 
density. Since Wg + Wi > 0, Qx diverges for large z. Also the BBN constraint (A) is not 
satisfied. 

III. Padmanabhan Parametrization(Jassal, Bagla and Padmanabhan,2004), 

WiZ 

For z » l,w{z) = Wo- Also ^ = ^^jfj^^- Therefore 'w{z) increases from 2; = to 2; = 1 
thereafter starts decreasing. Eq.(3) gives 

Qx = 0^(l + ^)^(^+-)xexp[^^^] (10) 

Best fit values for SNe 'gold set' Wo — —1.9,1^1 = 6.6, — 0.3, 2;^ = 0.3 give Qxi^r) = 
0.583 whereas constraint (B) demands Qx{zt) — 0.457 i.e. the parametrization gives a 
+27%deviation from the acceptable value of dark energy density. The BBN constraint (B) 
is satisfied since for z » 1, ^x{z) ~ (1 + 2;)''^(i+'"°) and 1 + Wo < 0. 

Therefore, Padmanabhan parametrization is compatible with astrophysical constraints 
and integrated tracking. 
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4 Hannestad Parametrization 

Let us now consider Hannestad parametrization (Hannestad and Mortsell,2004) which in- 
volves 4-parameters. 

W[Z) — WoWi = — ; 1 / -I , — (11) 

where a — {1 + z)~^, ar is the scale factor at the transition redshift zt, wq and wi are 
the asymptotic values of ^(2;) as 1 + z — > and z — > 00 respectively; n determines the 
rapidity of transition. 

It follows that the equation of state at transition epoch is given by 

wiaT) = (12) 

and the equation of state at the present epoch [z — 0) is given by 

Out of the 4 parameters, Wo and w\ are taken to represent the asymptotic values and w* 
and n are left free; alternatively wi and w* may be given pre-assigned values and Wo and n 
taken as free parameters to match with the observations. 
Differentiating Eq.(ll) with respect to z, we get 

_^dw {wi-Wo)n{l + zY-^d^ 

W — (14) 

dz [l + {aT/aY][wx + Wo{aT/aY] ^ ' 



Hence the equation of state of dark energy varies with the redshift unless W\ 

dz 



w{z) takes asymptotic values l + z—^Qoil + z^oo. The variation gradient |^| varies 
directly as \w\^ \wi — Wo\^n and inversely as aj.(l + z)^. 
At the present epoch(2; = 0) 

dw WoWi{wi — Wo)na^ , . 

dz {wi + Wottj^y 

At transition epoch(2; = zt) 



dw WoWi{wi — Wo)naJ}, 
dz '^iwi + Wo) 



(16) 



Note that changes sign in going from z = to z = zt ii Wo < 0,wi < and Wi — Wo > 0. 
Inserting for w{z) from Eq.(ll) into Eq.(3), we get 
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For z » 1 



Taking the best fit values for tlie 'gold set' SNe (Hannestad,2004), 

wi — —OA; Wo — — 1.8;n = 3.41 with a prior = 0.38, we have w{zt) — —0.654, 
by Eq.(12). 

The fourth parameter w* is still free; we can choose w* = —0.95, a value favored 
by observations, then Eq.(13) yields = 0.367. Substituting for Zt in Eq.(17), we 
find Qx{zt) = 0.68 whereas constraint (B) demands Qx{zt) = ~ 3w{zt) ~ 0.509. Thus 
Hannestad parametrization gives +33.3% deviation from the acceptable value of dark energy 
density. Again, it is obvious from Eq.(18) that the dark energy density parameter diverges 
for large values of redshift since 1 + Wi > for the best fit values given above; as such 
constraint (B) is not satisfied. 

5 Interpolation of w(z) for Slowly Varying Equation of 
State 

If we go by the analysis of Riess et al(2004) and Jassal ct al(2004), the current SNe data 
is inconsistent with rapid evolution of dark energy. Therefore, we can use 'integrating 
tracking' and interpolation techniques applicable to slow time varying equation of state 
(Johri,Pramana,2002), It was shown (Johri, Class. & Quantum Gravity,2002) that the scalar 
fields with slowly varying equation of state which satisfy tracking criteria, arc compatible 
with astrophysical constraints (A) and (B) outlined under 'integrated tracking'. According 
to our analysis (Johri, Pramana,2002), based on integrated tracking, we have 

Qx < 0.14, w ~ -0.035 at z = 10^° 
Ox = 0.66, w = -0.77 at z = 
Vlx = 0.5, w = -0.66 at zt = 0.414 

The above data is consistent with the transition redshift given by Riess et al(2004) and 
{^Ix)bbn given by Cyburt et al(2004) To apply Hannestad parametrization, if we choose 
asymptotic values wi = — 0.035, = —1 (assuming quintessence dark energy), we get 
w{zt) = —0.676 which yields zt — 0.405 and ilxizr) — 0.493. These values are remarkably 
close to the interpolated data given above. 

6 Conclusion 

We have investigated the behavior of evolving dark energy by taking parametric represen- 
tation of the equation of state. Out of the 2-parameter models of the equation of state, 
Padmanabhan parametrization is found to be valid and compatible with the astrophysical 
constraints over a wide range of redshift. 
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Bassett et al (2004) have discussed the hmitations of 2-index parametrizations. Most of 
them track well at low rcdshifts z < 0.2 but if we explore dark energy for regions z > 1 
and w << —1 (beyond quintessence models), very rapid variation in w{z) can be found. 
Higher order parametrizations are more suitable for probing the nature of dark energy and 
its evolution since more parameters give more freedom to fit in observational data but at the 
same time it gives rise to more degeneracies in the determination of the parameters, 

Hannestad 4-parameter model of the equation of state provides a well- behaved repre- 
sentation of dark energy evolution over a long range of redshift but it fails to satisfy the 
BBN constraint for dark energy density for quintessence models «; -|- 1 > at nucleosyn- 
thesis epoch around z — 10^°. However, if the weak energy condition (WEC) is relaxed to 
admit phantom models {w < —1), which are more favored by SNe and galaxy clustering 
observations(Caldwell et al,2003), not only BBN constraint is satisfied but as an additional 
bonus, Hannestad parametrization admits rapid variation of the equation of state as well. 
Since there is no lower bound on w{z) for the phantom models, flx{z) also decreases steeply 
in the early universe if 1 -|- is negative in Eq.(18). 
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